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DPOAE MEASUREMENT SYSTEM AND 

METHOD 

Abstract 

A new method of DPOAE signal level measurement based on a recently developed family 
of nonlinear adaptive algorithms is presented. Three core units are employed to extract 
and measure the two artifacts and the DPOAE signal. Each core unit, which can be any 
of the three members of a new family of signal processing algorithms, has the capability 
to lock on and track the time variations of a specified sinusoidal component of its input 
signal The core units have very simple structure while maintaining very good perfor- 
mance in noisy environments in terms of stability and convergence speed. Performance of 
the present method is demonstrated with the aid of computer simulations. The present 
method features structural simplicity which renders it suitable for both software and 
hardware implementation. It offers extremely high noise immunity and requires shorter 
measurement time than conventional methods. 
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Introduction 

Distortion product otoacoustic emissions (DPOAEs) are very low level stimulated acoustic 
responses to two pure tones presented to the ear canal. DPOAE measurement provides 
an objective non-invasive measure of peripheral auditory function and is used for hearing 
assessment [1]. DPOAE screening is becoming a standard clinical practice to predict 
potential sensorineural hearing loss especially in newborns. 

DPOAEs have been recognized for a number years [2, 3]. However, DPOAE measurement 
is considered an active area of research because of the challenging nature of the signal 
processing task. To address the ever-increasing demand for high performance DPOAE 
measurement methods, a number of signal processing algorithms have been present in 
recent years [1, 4-6]. With the availability of the powerful computational tools such as 
digital signal processors (DSPs), commercial medical equipment dedicated to DPOAE 
measurement are becoming available [7-14]. 

In this type of otoacoustic test, two pure tones with frequencies /i and /2 are presented 
to the cochlea. For best results, / 2 is usually chosen as 1.2 f\. Due to the non-linearity 
of ear, a very low level of distortion product of frequency 2/i - / 2 is generated in normal 
ears. The level of such DPOAE signal is a measure of functionality of the ear. Estimation 
of such a weak signal buried under two strong artifacts in a potentially noisy background 
is a challenging signal processing problem. 

Conventionally, fast Fourier transform (FFT) is used as the main signal processing tool to 
estimate the level of DPOAE signals. Application of FFT in this problem has a number 
of shortcomings among which long measurement time is most pronounced. Such long 
measurement time is usually required for acquisition of more data which, when averaged, 
reduce the overall background noise effect. Unreliability of the measurements is another 
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problem of FFT-based methods and is a direct effect of the sensitivity of the FFT-based 
methods to the background noise. In addition to the need to increase the measurement 
time, the tests are usually required to be conducted in low noise environments such as 
sound-proof booths. 

In attempt to devise high performance DPOAE estimation techniques, linear adaptive 
signal processing techniques have been employed [1, 4, 14]. Such techniques generally 
offer better performance in terms of measurement time which may be interpreted as 
higher noise immunity of adaptive techniques compared to FFT. However, the need for 
sound-proof examination rooms is not obviated with such techniques. 

This document presents a method of measurement of DPOAE signals based on a new 
signal processing technique in which individual sinusoidal components of a given signal 
are extracted and their variations are adaptively tracked over time. The core units are 
presented in a separate disclosure [15]. The present DPOAE estimation method consists 
of three units which can be any of the three such algorithms. The two artifacts are first 
extracted by two units and are subtracted from the input signal to generate an input 
signal of which DPOAE signal has a higher relative portion. Such a signal is then fed 
to another core unit which estimates the level of DPOAE signal. Superior performance 
of the present technique in terms of noise immunity and fast measurements is illustrated 
with the aid of computer simulations. 

In order to provide a brief background on the subject, a general review of the structure 
of a generic DPOAE measurement device is presented in the next section. It will be seen 
that the heart of such an apparatus is the signal processing subsystem which, in the final 
analysis, determines the performance of the overall system. Section presents the present 
scheme for such signal processing module. The core units which are the building blocks 
of the present signal processing scheme are introduced in section . Performance of the 
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present method is demonstrated in section . 

Structure of a Generic DPOAE Measurement Device 

In this section a brief overview of the structure of a typical DPOAE detection system is 
provided. Figure 1 shows the generic block diagram of a DPOAE measurement device. 
It consists of three main modules: the data acquisition/transducers module, the signal 
processing module and the display. 

Data acquisition unit is the medium between the processing unit and the probe which 
transmits and receives acoustic signals in the audio range. Components of the compound 
data acquisition/transducers module are illustrated in more detail in Figure 2. One of the 
main functions of this module is to convert digital signals produced by the signal process- 
ing module to analog signals which are then conditioned and converted to audio signals. 
Conditioning of the signals in this case may or may not include filtering. Conversely, the 
audio signals recorded by the probe are conditioned and converted to digital signals to be 
processed by the signal processing module. 

The signal processing module is the heart of the system which produces the digital form 
of the artifacts and extracts and measures the DPOAE signal. A DSP, or if the com- 
putational/architectural demand is low even a microcontroller, can be employed as the 
hardware platform of this unit. Signal processing is embedded as the software in such a 
hardware platform. Alternatively, and provided that the complexity of the signal process- 
ing algorithms remains low, signal processing unit may be implemented solely in hardware 
using programmable logic array (PLA) or field programmable gate array (FPGA) tech- 
nology. In an ideal case, namely when the signal processing algorithm is not excessively 
complex, the hardware does not require a PC for its operation; however, interfacing to a 
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PC is usually provisioned for data management. 

The display unit is the interface between the device and the operator. It can be a simple 
LED/LCD and/or a small printer. 

Present Technique 

Figure 3 shows the main functions of the software embedded in the signal processing 
module. The software is essentially responsible for the generation of the artifact signals 
and extraction of DPOAE as well as management of input/output data. As discussed 
before, the significance of the present work is in the introduction of a signal processing 
technique for the extraction and measurement of DPOAE signals. 

The present signal processing scheme employs three core units to construct a high perfor- 
mance DPOAE extraction module. The three possible structures of such building block 
units are presented in section . Each core unit is capable of focusing on and extracting 
a pre-specified sinusoidal component of its input signal which may contain many other 
components including noise. More importantly, they can effectively follow variations in 
the amplitude, phase (and frequency) of the extracted sinusoidal component. Although 
the underlying mathematics ensuring stability and performance of such core units is very 
complex, the structure of the core units remain extremely simple. They are found to be 
very robust with respect to variations in the internal settings as well as external conditions 
and exhibit superior performance over existing linear adaptive and FFT-based algorithms. 

The input signal is often assumed to consist of two pure sinusoids with frequencies f\ and 
f 2 at a very high level (usually about 60-70 dB) and a very low level DPOAE 2/ x - f 2 at 
about 0 dB. It is polluted by a noise usually considered to be at about -10 dB level. The 
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noise in fact represents the totality of all undesirable signals that may be present in the 
environment in which the examination is being conducted as well as unavoidable white 
Gaussian noise. Because of excessive degree of pollution (artifacts and noise), one single 
core unit set to extract the DPOAE signal out of the input exhibits poor performance. 
Different arrangements were studied to construct a high performance architecture. The 
most successful configuration is shown in Figure 4. Three core units are employed. The 
first two core units are set to extract the artifacts. They effectively do so with very small 
errors. The extracted artifacts are then subtracted from the input to produce a signal, of 
which DPOAE has a higher relative portion. The third core unit is then set to extract 
DPOAE. 

Review of the Core Unit 

This section reviews the mathematical structure and properties of the core units which are 
the basis for the present DPOAE measurement method of section . A complete description 
of the signal processing core is provided in [15]. Let u(t) denote a signal comprising a 
number of individual sinusoidal components and noise, expressed by 

N 

t*(*) = £4fe sin U) 

where fa = o^t + 0* is the total phase, and n(t) denotes the total noise imposed on the 
signal. The objective is to find a scheme for estimating a certain component of such input 
signal as fast and accurate as possible; a scheme which should not be sensitive to the noise 
and the potential time variations of the input signal. Simplicity of the structure, for the 
sake of practical feasibility, is desirable. 

Let M. be a manifold containing all pure sinusoidal signals defined as 
M = {y(t, 0) = Ox sin(0 2 t + 0 3 )| ft G 
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where 0 = [0i,02,*03] T is the vector of parameters which belongs to the parameters space 

and T denotes matrix transposition. To extract a certain sinusoidal component of u(t), the 
solution has to be an orthogonal projection of u(t) onto the manifold M, or equivalently 
has to be an optimum 0 which minimizes a distance function d between y(t, 0) and u(t), 
i.e., 

0 opt = aigmmd(y(t,0),u(t)). 
In least squares method, d is the instantaneous distance function 

d\t,0)^{u{t)-y{t,0)f^e{tf. 

The parameter vector 0 is estimated using the gradient descent method, 

where the positive diagonal matrix \i is the algorithm regulating constant. It controls the 
convergence rate as well as the stability of the algorithm. 

Following the steps outlined above, three possible sets of differential equations are obtained 
depending on the definition of the form of the parameter 0. Complete derivations are 
presented in [15]. The governing sets of equations of all three members of this family of 
the algorithms can be written under the same formulation as 

A = //iesin<£, (2) 

lj = /^eAcos^, (3) 

4> = F(A 5 u,0,e), (4) 

y{t) = Asin<£, (5) 

e(t) = u(t)-y{t\ (6) 

where y(t) is the output and the function F(A y u, 0, e) can take any of the following forms: 
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• Type 1: F{A,u,<f> y e) = ii 2 e A cos <j> + u 0 , 

• Type 2: jF(A,u;,<£,e) =tu, 

• Type 3: F[A, cj, 0, e) = p^eA cos + uj 

where u; 0 is a fixed center frequency. The dynamics resulting from the first of these forms 
can be expressed in a more concise way as 

A = /ixesin^, 

<j> = foe A cos <j> + u; 0> 
y(t) = Asin^, 
e(t) = u(t)-y(t). 

It has been shown that each of the three dynamical systems represented by the above sets 
of differential equations possesses a unique asymptotically stable periodic orbit which lies 
in a neighborhood of the orbit associated with the desired component of the function u(t). 
In terms of the engineering performance of the system, this indicates that the output of 
the system y(t) = Asin0 will approach a sinusoidal component of the input signal u(t). 
Moreover, the slow variations of parameters in u(t) are tolerated by the system. 

Figure 5 shows implementation of the algorithms in the form of composition of simple 
blocks suitable for schematic software development tools. Numerically, a possible way of 
writing the set of equations governing the present algorithm in discrete form, which can 
be readily used in any programming language, is 

A[n -hi] = A[n) -h T 5 /zie[n] sin <j>[n] , 
w[n -hi] = u)[n) + T s /i 2 e[n]i4[n] cos 0[n], 
#x+l] = 0[n] + T a F(A[n],u;[n],0[n],e[n]), 
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y[n) = A[n]sin<£[n], 
e[n] = u[n]-y{n] 

where a first order approximation for derivatives is assumed, T 3 is the sampling time and 
n is the index of iteration. 

The dynamics of each of the three algorithms presents a notch filter in the sense that it 
extracts (i.e. lets pass through) one specific sinusoidal component and rejects all other 
components including noise. It is adaptive in the sense that the notch filter accommodates 
variations of the characteristics of the desired output over time. The center frequency of 
such adaptive notch filter is specified by the value of u) 0 in the first type or the initial 
condition of the variable iv in the second and third type. 

A physical system approaches its steady state in about a few times its time constant 
r = jj where [i represents the eigenvalue of the system. The eigenvalues of the presented 
dynamical systems are — ^ and — *?M 2 , which determine the speed of the system in 
approaching the solution. For example, a set of values of ^1=^2 = 10000 should cause 
the introduced dynamical systems to reach their steady state in about 5r = 1 ms for an 
input signal of unit amplitude. 

It has been observed that all three types of the algorithms introduced in this section can 
be employed as the core units of the present DPOAE measurement scheme of Figure 4. 
For the purpose of demonstration of the performance of the algorithms and the DPOAE 
measurement scheme, type 1 algorithm is used throughout this document as the building 
block. Figure 6 shows the performance of one such unit in the extraction of a specified 
sinusoidal component of its input signal. As noted before, the specification of the desired 
sinusoidal component is achieved by means of appropriate setting of parameter u 0 = 2nf 0 
which determines the frequency of the component of interest. The input signal is a clean 
sinusoid of frequency / = 3000 Hz with a random constant phase. All the initial conditions 
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are taken as zero. As would be expected, the convergence is achieved within about 1 ms. 

When the input signal consists of a single pure sinusoid which is to be extracted, the 
values of the parameters fi\ and fjL 2 can be chosen very large so as to increase the speed of 
convergence without any trade-off with accuracy. However, as the degree of the pollution 
in the input signal, which may be quantitatively represented by total harmonic distortion 
(THD) or signal to noise ratio (SNR) increases, there will always be a trade-off between 
the speed and accuracy. Figure 7 shows the performance of the single , unit set as before 
for a sinusoidal input polluted by a white Gaussian noise with SNR=20 dB. With this 
parameter setting, a noise of 100% magnitude in the input introduces about 5% error in 
the estimated value of the amplitude. The same experiment is now repeated with a much 
smaller values of parameters \i\ and ^ Figure 8 shows the result of this case. It is clear 
that a reduction of parameters by a factor of 60 results in a convergence time of about 60 
times longer. For this, one gains a reduction in the incurred error of a factor of about 10. 

Another factor to be considered in the assignment of values to parameters /ii and /i 2 is 
the amplitude of the input signal. As noted before, the eigenvalues of the algorithm are 
— ^ and -*fA 2 . The first one is merrily independent of the input signal, but the second 
one is dependant on the amplitude of the input signal. In order to achieve the same speed 
of convergence when the amplitude of the input signal is reduced, one has to increase 
the value of the parameter fi 2 by the square of the reduction factor. Figure 9 shows the 
performance of the algorithm when this provision is made which verifies the theoretical 
prediction. 
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Performance 

This section presents the performance of the present technique by the use of computer 
simulations. Three core units of the first type are used in the arrangement of Figure 4 to 
construct the present DPOAE measurement scheme. All the initial conditions are taken to 
be zero. The first two units which are employed to extract the artifacts must do so with as 
small an error as possible to produce an input to the third unit containing DPOAE signal 
as a main component, otherwise the third unit will not be able to extract DPOAE signal 
with sufficient accuracy. Therefore, the parameters are chosen relatively small at \i\ = 100 
and fa = 50. With these values, it is expected that the convergence is achieved within 
about 100 ras. The same setting was also found suitable for the third unit. Two forces 
oppose each other in the selection criteria for the values of the parameters of the third 
unit. On the one hand, since SNR is very low for the third unit (DPOAE being taken as 
the signal and the total residual errors of the first two units as the noise), the parameters 
must be very small. On the other hand, since the input has a very low level, the value of 
/x 2 must be large. These two forces seem to cancel each other , for the case of \ii so that 
the original value of 50 remains a good choice. As for the value of one expects that it 
should be small to achieve acceptable accuracy in the extracted signal. This is true, but 
since only the amplitude of the extracted signal (i.e. the level of DPOAE) is of interest, 
one can substantially increase The accuracy in the estimated extracted signal (and 
not necessarily its amplitude) is thus traded off which in this case is not of importance. 
The estimation of the amplitude remains accurate while the convergence speed is high 
enough. At any rate, a rough setting of the values of parameters /ii and /x 2 axe to be done. 
The core algorithms have been shown to be very insensitive with regard to the settings of 
the parameters; variations as much as 50% in the values of these parameters seem to have 
practically no effect on the performance. The value of the amplitude directly estimated 
by the third unit was used as the estimate of the level of DPOAE. 
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The input signal is synthesized by the summation of two artifacts at frequencies f\ = 
3000 Hz and f 2 = 1.2/i = 3600 Hz and the DPOAE signal at 2f x - f 2 = 2400 Hz. The 
amplitude of the artifacts are taken as 1 v and the DPOAE level is 0.5 mv. With these 
values, the level of the artifacts is set at about 66 dB higher than that of DPOAE. Figure 
10 shows the performance of the present technique. It is confirmed that the measurement 
is completed in about 100 ms. 

To present the dynamic performance of the present algorithm, the response of the algo- 
rithm to a step change in the level of DPOAE signal is presented in Figure 11. It is clear 
that the algorithm follows the variations in the level of DPOAE at a sufficiently high 
speed. 

A number of measurements can be made consecutively. Figure 12 shows such an example. 
In each case, the frequency of the extracted DPOAE signal is shown. As the frequency 
decreases, less number of cycles are available in a given time interval. For example, there 
are 160 cycles available in 100 ms of a signal of frequency / = 1600 Hz whereas in the 
case of a signal of frequency / = 4000 Hz the available number of cycles are 400 in the 
same time interval. This explains why the estimation of DPOAE level is better achieved 
at higher frequencies. It should be noted that within the range of frequencies of interest 
for DPOAEs, the algorithm provides acceptable estimates in about 100 — 200 ms. 

As noted before, one of the main features of the employed core algorithms is their noise 
immunity. Therefore, one expects that the overall DPOAE estimation scheme be immune 
to background noise. Usually, otoacoustic tests are conducted in quiet examination rooms 
so that the level of noise floor is low enough for good results. Under such controlled 
conditions, a noise floor of about 10 dB lower than DPOAE signal (SNR=10 dB) exists. 
It is obvious that an algorithm which is able to extract a DPOAE signal highly immersed in 
background noise is desirable in that it alleviates the need for special quiet examination 
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facilities such as sound-proof booths. Figure 13 shows the performance of the present 
algorithm in a highly noisy environment where the level of the DPOAE signal is set to be 
at the level of the noise floor (SNR=0 dB). The incurred error is negligible. In Figure 14, 
the level of noise is elevated to be 10 times larger than the that of DPOAE signal itself. 
The incurred estimation error is less than ±15% which for most practical application is 
sufficiently good. The excellent noise immunity of the present algorithm not only obviates 
the need for sound-proof examination rooms, but also provides the way to reduce the level 
of the artifacts for more patient-friendly tests. 
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Figure 1: Generic block diagram of a DPOAE measurement device. 
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Figure 2: A more detailed diagram of data acquisition and probe units. 
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Figure 3: Main functions of the embedded softwaxe. 
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Figure 4: Block diagram of the signal processing scheme. 
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Figure 5: Block diagram implementation of the three possible core units of type 1 — 3. 
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Figure 6: Performance of a single unit in the extraction of a sinusoidal component and its 
amplitude. The dashed line waveform is the input. The parameters are set at = 12000 
and \i 2 = 6000. 



4 



CA 02346243 2001-05-04 




1.4 1.6 

Time (ms) 



Figure 7: Performance in the presence of noise. The parameters are set at fix = 12000 
and fa = 6000. The bottom figure is a magnified portion of the top figure. 
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Figure *8: Performance in the presence of noise. The parameters are set at = 200 and 
& = 100. The bottom figure is a magnified portion of the top figure. 
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Figure 9: Performance of a single unit in the extraction of a sinusoidal component and its 
amplitude when the amplitude of the input signal is reduced by factor of 10. The dashed 
line waveform is the input. The parameters are set at /zi = 12000 and ^ = 6000 x (10) 2 . 
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Figure 10: Performance of the DPOAE estimation method in measurement of a signal of 
0.5 mv amplitude. 
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Figure 11: Dynamic performance of the present method. The level of DPOAE undergoes 
a step change at t = 200 ms. 
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Figure 12: Three consecutive measurements with different frequencies. 



10 



CA 02346243 2001-05-04 



0.7 



0.6 




£0.48- 



UJ0.47 

0.46- 

°* 4 ?00 120 140 160 180 200 

Time (ms) 

Figure 13: Performance of the present method in estimation of DPOAE in a highly noisy 
background (SNR=0 dB). The bottom figure is a magnified portion of the top figure. 
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Figure 14: Performance of the present method in estimation of DPOAE in an extremely 
noisy background (SNR=-20 dB). The bottom figure is a magnified portion of the top 
figure. 
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